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We constrain a primordial magnetic field (PMF) generated during a phase transitions 
(PT) using big bang nucleosynthesis (BBN) bounds on the relativistic energy density. The 
amplitude of the PMF at large scales is determined by the spectral shape of the PMF 
spectrum outside its maximal correlation length scale. Even if the amplitude of the PMF 
at 1 Mpc is small, PT generated PMFs can leave observable signatures in the potentially 
detectable relic gravitational wave background if a large enough fraction (1 — 10%) of the 
thermal energy is converted into the PMF. 

1. Introduction 

A cosmological seed PMF (generated during or prior to the radiation-dominated epoch) 
has been proposed to explain the existence of observed ~ 10~ 6 — 10~ 5 Gauss (G) magnetic 
fields in galaxies and clusters [HI2]- To preserve approximate spatial isotropy a PMF has 
to be small and hence can be treated as a first order term in perturbation theory. In the 
standard cosmological model [E] the energy density parameter of a PMF, Q B = ps/Pcr, is 
significantly less than unity. Also, a PMF must be smaller than those observed in galaxies 
(10 -5 G), so f^/iQ < 10~ 4 where ho is the Hubble constant in units of 100 km s _1 Mpc -1 . 
Since the PMF energy density contributes to the radiation field, the BBN bound implies 
^Bh.Q < 2.4 x 10~ 6 . The ratio of pB and the energy density of radiation p ra( j is constant 
during cosmological evolution, if the PMF is not damped by an MHD (or other) process 
and so stays frozen to the plasma. Direct measurement of a cosmological MF is based on 
the Faraday rotation effect. A potential extension of this method, based on the rotation of 
the cosmic microwave background (CMB) polarization plane, appears promising [HIE]. I n 
addition, a PMF leaves imprints on the CMB temperature and polarization anisotropies 
(for a review see Ref . [ [6] ) . 



* tinatin@phys.ksu.edu 
t aleko@tevza.org 
•tratra@phys.ksu.edu 



2 



Kahniashvili et al. 



In this paper we consider cosmological PMFs generated by causal processes during 
phase transitions (PTs) such as the electroweak (EW) and QCD PTs [ [7]. The main 
parameters of interest are the temperature T* and the number of relativistic degrees of 
freedom g+ when the PMF is generated. We only use fundamental physical laws, such as 
conservation of energy and how the magnetic field interacts with the cosmological plasma 
through MHD turbulence, and do not make any assumption about the physical process 
leading to PMF generation. We discuss cosmological signatures of such a PMF, including 
the effects on the CMB temperature and polarization anisotropies and the production of 
gravitational waves (GWs). We employ natural units with h = 1 = c and gaussian units 
for electromagnetic quantities. 

2. Magnetic Field Spectrum 

The maximal correlation length Z max at generation of a causally generated PMF cannot 
exceed the Hubble radius if" 1 then, so 7 = Z max / ' H~ l < 1, where 7 can be associated 
with the number of PMF bubbles within the Hubble radius, N oc 7 s . The comoving 
length (measured today) corresponding to the Hubble radius is inversely proportional to 
the temperature T*, 



and is equal to 0.5 pc for the QCDPT (with g+ = 15 and T* = 0.15 Gev) and 6 x 10 4 pc 
for the EWPT (with g+ = 100 and T* = 100 GeV), and the comoving PMF correlation 
length £ max < X H . 

If generated prior to BBN, the maximal value of the PMF energy density must satisfy 
the BBN bound, i.e. the total energy density of the PMF at nucleosynthesis, pb(on), 
should not exceed 10% of the radiation energy density then, /w^n)- Since the ratio 
Ae/Prad is constant, the maximal comoving value of the effective PMF = y/8npB = 
8.4 x 10" 7 (100/V) 1/6 G, if no PMF damping occurs before BBN. Even if the PMF energy 
is converted to another field contributing to the radiation (for example, GWs [M, E]), there 
is only pe(a*) magnetic energy available. The next issue is to determine how this energy 
is distributed at different wavelengths, and the comoving PMF at a given comoving length 
scale A. Of course, for a scale-invariant [[TO] or homogeneous PMF the limit remains the 
same at any scale. Note that the maximal value of the PMF (from the BBN bound) is 
independent of the temperature at generation T*, and depends only very weakly on the 
number of relativistic degrees of freedom then. 

Pb can be viewed as the magnetic energy density injected into the cosmological plasma 
at the comoving length scale Ao, the size of the largest magnetic eddy. After generation, 
PMF evolution (during the PT) depends crucially on the length scale under consideration. 
If the relevant time scales are shorter than H~ l we can neglect the expansion of the 
Universe. In this case, we must distinguish three sub-Hubble-radius regimes: ku < k < k 
(where ku = 2tt/\h and k = 27i/\ ; the large scale decay regime), ko < k < kn with 
kr> = 2n/\ D the damping wavenumber scale related to plasma properties (the turbulence 
regime), and k > kr> (the viscous damping regime). The interaction of the PMF with the 
plasma, and as a consequence the dynamics of the PMF, is sensitive to the presence of 
magnetic helicity (see Refs. [[11] for magnetic helicity generation mechanisms). 
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The magnetic energy EM{k, t) and helicity HM{k, t) density power spectra are related to 
the magnetic energy and helicity densities through Euif) = J^°dkE M (k,t) and Tiuif) = 
J^°dkH M (k,t). The magnetic correlation length — [J™dkk~ l EM{k, t)]/£M(t) cor- 

responds to the largest eddy length scale. All configurations of the MF must satisfy 
the "realizability condition" [[12]: |Hm(^)| < ^ M {t)£ M {t). Also, the velocity energy 
density spectrum Ex(k,t) is related to the kinetic energy of turbulent motions through 



2.1. Magnetic Field Spatial Structure 

We first consider the non-helical case. For large enough Reynolds number the magnetic 
energy is re-distributed by a Kolmogoroff turbulence direct cascade. From the analogy 
between the Kolmogoroff laws for hydrodynamic and magnetic turbulence, the magnetic 
energy dissipation comoving rate per unit enthalpy is Em — (2/3) 3 ^ 2 kov\, with va = 
a/1.5pb/ p ra d being the effective Alfven velocity corresponding to the total fluid- injected 
PMF energy, i.e. Em = fco(PB/Prad) 3 ^ 2 - In the absence of PMF damping at k = k , 
for Kolmogoroff turbulence, -Em(^) = CmPb&o" 1 ^ ^ 3 when k < k < ko-, where Cm is 
a constant of order unity and k = k/k®. At large scales when k < k$ we model the 
PMF energy spectrum by a power law, E(k) oc k a . Requiring continuity of the PMF 
spectrum at k = ko, E{k) = CMPsk^k 06 for k < k$. It is natural to assume that the MF 
energy injection scale Ao is the same as the maximal correlation length of the PMF, i.e. 

Ao — ^max fl o/ <V 

The largest scale MF energy density spectral index a has been much discussed. Hogan [ 
[13] requires causality of the field and argues that the PMF energy density spectrum must 
be white noise for scales larger than the causal horizon, Ao- This corresponds to a = 2. 
Durrer and Caprini [[H] claim that this violates the divergence-free PMF requirement and 
instead demand a = 4. Both of these spectra, a = 2 (Saffman) and a = 4 (Batchelor)Q 
are well known in the turbulence literature and, as discussed in Ref. [US], their realization 
depends on initial conditions. Another possibility is Kazantsev's a = 3/2 value, which 
can be rapidly achieved during the turbulence decay process discussed in Refs. [US]. To 
keep the analysis as general as possible, we keep a arbitrary as much as possible. 

Requiring E M < p B we obtain Cm < 1.5 (a + l)/(3a + 5). With p B < 0.1p ra( j and ne- 
glecting the contribution to the energy density from scales smaller than the damping scale 
Xd, for the maximal allowed value of Cm, we have CmPb < 2.81 x lCT^CA/.maxllOO/g*) 1 / 3 
G 2 , and the magnetic energy spectrum 



Note the damping wavenumber kn is determined by the Reynolds number Re ^> 1 of the 
fluid during the PT, as k D = k Re 3 ^. 

Defining B\ as a smoothed PMF over a sphere of radius A (A > Ao) we have for the MF 
energy density on scales larger than the maximal correlation length, £j^ s = dkE M (k) 
= Bl(k Q X) a+1 /[8nT(0.5a + 1.5)] where V is the Euler Gamma function, [H?]. This leads 
to the upper bound B\ oc X~^ a+1 ^ 2 shown in Fig. 1. The PMF limits shown in Fig. 1 
does not account for the evolution and damping of the PMF during the expansion of the 

A E r~j k is sometimes called the von Karman spectrum. 



£ K (t) = f™dkE K (k,t). 




max 




(2) 
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Figure 1. The maximal allowed value B\ for a PMF generated during the EWPT with 
T* = 100 Gev, </* = 100, and 7 = 0.01. 
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Universe, addressing only the spatial structure of the PMF at large scales constrained 
by the BBN bound on the relativistic component energy density. Based on dimensional 
analysis we address the time evolution of the PMF in the next subsection, with more 
precise computations presented in Ref. [US]. 

2.2. Magnetic Field Temporal Characteristics 

The PMF spectrum is characterized not only by its spatial distribution, but also by 
its characteristic times: i) the largest size eddy turn-over time r ~ Iq/va, which can 
also be used to determine the minimal duration of the source needed to justify use of 
the stationary turbulence approximation [[THjJfl and, ii) the turbulence cascade time-scale 
Td c and the large-scale turbulence decay time t\ s . If the source duration time is short 
compared to the Hubble time H~ l we can neglect the expansion of the Universe. For the 
EWPT (7 ~ 0.01), and for reasonable moderate values of the Alfven velocity (va < 0.3), 
the maximal size magnetic eddy turn-over time scale To — , jH~ 1 /va is significantly shorter 
than the Hubble time H~ l . Also, t& c is determined by the dissipation rate Em- To proceed 
we specify the time decorrelation function f(r](k),r), as fdc[v(k),T~] = exp[— nrf {k)T 2 / A\ 
for ko < k < kr, [ 121] . where r](k) = e l ^k 2 ^ / ' . For non-helical turbulence, by 
considering the largest size magnetic eddy decorrelation, r^c — 0.5r , thus the direct 
cascade time scale is much shorter than the Hubble time, and the assumption made 
above to neglect the energy density for A < A is justified, as is the assumption to neglect 
the expansion of the Universe during the direct cascade. 

The other characteristic time is related to the decay of large-scale turbulence. Specific to 
this process is that there is no magnetic or hydrodynamic turbulence production source 
and free decay occurs. In this case here we may adopt the grid turbulence decay law 
k(t) oc t~ 1,3 (see e.g. [E2])- This can be motivated by noting that the correlation length 
and the Hubble radius set natural length scales which are the analogue of the grid size 
in laboratory turbulence. Assuming that k is time independent during the decay and 
£m(£) ^ ^~ 2 ) we get E M (k,t) oc k a t~ - 7 when k < 1. Figure 2 shows the spatial and 
temporal surface of freely decaying turbulence at large scales ignoring MHD dynamo 
effects. The time scale is normalized to the largest size magnetic eddy turn-over time r . 
For scales k in the range k Q < k < kjj, temporal decorrelation occurs in a time interval 
much shorter than the time scale related to the free decay of turbulence. Thus the free 
decay timescale is important only at large scales. The grid turbulence analogy implies 
that the PMF spectral energy density is decreased by several orders of magnitudes within 
a period of IOto (well within the EWPT Hubble time H~ l ), indicating that the expansion 
of the Universe can be neglected. After free decay reduces the MF power, the PMF on 
scales k <C k can be treated as a MF that is unaffected by turbulence. Of course, when 
considering realistic cosmological turbulence, in contrast to the laboratory turbulence 
case, the evolution of the fluid viscosity must be taken into account, which can change 
the PMF correlation length and energy density scaling laws shown here, see also Ref. [ 
[23] for a different model of free-decaying MHD turbulence. We also argue [ ITS] that free 

5 Since, for the developed turbulence case, the magnetic energy proper dissipation rate per unit enthalpy, 
£m = CM^o/a-*, must be approximately equal to the mean energy input rate per unit enthalpy of the 
source, i.e. the turbulence cascade time scale Td c — 27r£^ 1 (pB/w ra( j) — tq where w ra( j is the radiation 
enthalpy [|20]. 
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Figure 2. Spectral energy surface of the large-scale non-helical decaying turbulence 
logEjtfik, t). Free decay of turbulence starts at t/ro = 1. Exponential decorrelation 
occurs in the inertial region where k > fco- In a lOro time period the spectral energy 
maximum drops by 5 orders of magnitude. 



decay laws are strongly initial condition dependent, i.e. the free decay law depends on how 
the PMF was generated [[?]; whether the PMF was generated through bubble collisions, 
leading to l ~ v^P' 1 with t> fc the bubble wall velocity and (3 the bubble nucleation rate 
parameter ((3 ~ lOOif* for the EWPT) [123], or if the PMF was present prior to the PT [ 
[TO] . Another uncertainty comes from the kinetic (vortical) energy density spectrum Ex [ 



The presence of even a small amount of magnetic helicity substantially affects PMF 
evolution [ [25], [26], [23 [28]. If there is only a little magnetic helicity, first a direct cascade 
develops. At the end of this first stage the turbulence relaxes to a maximally helical 
state [[2Hl[2Z] that satisfies < 2^(t)^(t) and the second inverse-cascade stage 

starts. Conservation of magnetic helicity implies that the magnetic energy density decays 
in inverse proportion to the correlation length growth during the inverse cascade. In 
contrast for the case of well-established non-helical turbulence, the effect of magnetic 
helicity is still under discussion [[212, [251 [23 I2H]- The main point of debate is related to 
the magnetic correlation length growth rate during the inverse cascade, i.e. £m(£) oc t n z, 
where the index is argued to be 1/2 [ [25], [2E] or 2/3 [[23 [21] • The total magnetic 
energy density oc t~ n ^, and the decay of large-scale kinetic energy Exit) (and as 

a consequence the ratio between the magnetic and kinetic energy densities Em{P) / 'Exit)) 
are sensitive to n^. In particular, Refs. [ [25] [26] argue that Exit) oc i -1 , implying a 
faster decay of kinetic energy at large scales, while the results of Refs. [ [23 I2H] lead to 
a constant Em($)/£k(P) within the inverse cascade and Exit) oc t -2 / 3 . In both these 
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models the magnetic turbulent energy density significantly decays on the phase transition 
timescale H~ x , see Fig. 3. 

3. Results and Discussion 

The limits on the PMF at large scales are much stronger for the non-helical turbulence 
case; without the decay the constraint at zero redshift on 1 Mpc is 10~ 28 G for a = 4 
for an EWPT generated PMF, see Fig. 1 (also see Ref. [ [23]). In the a = 3/2 case, 
the PMF can reach values of order 10~ 12 — 10~ n G that are required for seed MFs that 
might be able to explain observed MFs in galaxies and clusters [ [29] . Of course, large- 
scale decay of turbulence will strengthen these limits for both the non-helical and helical 
cases. On the other hand, accounting for large-scale PMF decay the BBN bound does not 
imply Pb < 0.1p ra d when the PMF is generated. However, there is another requirement: 
the PMF cannot be the only component during the radiation dominated epoch, thus 
Pb/ Prad(^*) < 1- Even though our analysis is preliminary, it seems that a PT generated 
PMF requires an effective amplification mechanism (such as a dynamo), or a specific 
initial condition, to act as a viable seed field for observed MFs in galaxies and clusters. 
We will address this issue in future work [ [18] . 

A PT generated PMF may have observable cosmological signatures. In particular a 
PMF induces CMB anisotropies. Usually when considering PMF limits from CMB data 
one refers to the amplitude of the smoothed PMF, B\, on large scales A ~ 1 Mpc [E]- 
On the other hand, when computing PMF induced CMB temperature and polarization 
anisotropy power spectra one finds C\ oc (B\\ a+1 ) 2 [ [XTl EH [31], and when considering 
Faraday rotation of the CMB polarization plane the rotation angle and the resulting B- 
polarization power spectra are oc B\\ a+1 [|5]. From the definition of B\ above, it is clear 
that PMF imprints on CMB fluctuations are determined by Vis (or Q 2 B ), also see Refs. [ 

E2JE3]. 

On the other hand, the MF energy density evolves because of cosmological expansion, 
resulting in the time dependence of the damping wavenumber kp. During a PT that 
generates the MF, or when the MF starts to interact with the plasma (if it had been 
generated prior to the PT), the total MF energy density at large scales is determined by 
the wavenumber of the peak of the MF spectral energy density, k . The value of this 
peak is fixed by the maximal correlation length of the MF at the phase transition, and it 
is independent of future evolution and dissipation processes. In other word, during PTs 
all wavenumbers k < k$ contribute to the total MF energy density. On the other hand, 
near the last scattering surface some modes have been damped and dissipated and thus 
the total energy density accounts for all wavenumbers k < k^. k^ is not determined 
by MF dissipation by viscosity on small scales during the PT, but by other processes 
that occur during cosmological expansion. Note that kr> < Refs. [ [30j [32] study 
the damping of a homogeneous MF assuming the main dissipation process is Alfven wave 
viscosity, resulting in k^ 1 ~ L$va, with Ls being the Silk damping scale at recombination 
and va determined by the background homogeneous PMF. For the case of a stochastic 
PMF, Ref. [EI] defines v A through the effective background MF, B eS = J B A (A/A D ) (a+1)/2 , 
smoothed over the length scale A > Ad, so ko becomes B\, A, and a dependent. Such 
a description appears inconsistent with the picture of Alfven wave induced dissipation, 
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Figure 3. Spectral energy surfaces of the large-scale helical turbulence \ogE M (k,t). The 
upper and lower panels correspond to the = 1/2 and = 2/3 decay laws, respectively. 
In a 100r time period the maximal spectral energy drops by 3 and 5 orders of magnitude 
in the two cases, respectively. 
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since the effective background field when a > is significantly smaller than that associated 
with the fluctuating Alfven wave MF. As a result the damping scale <C k Q and only a 
small part of the initial MF energy contributes to the CMB anisotropies. Accounting for 
this, a PMF with a = 4 and amplitude of order 10~ 13 G at 1 Mpc contributes to CMB 
fluctuations, while a scale-invariant PMF with a ~ — 1 and amplitude larger than 1CU 9 G 
might leave observable CMB anisotropy traces. On the other hand, ignoring Alfven wave 
induced damping (and not determining va through the amplitude of the PMF at scales 
corresponding to Lg) we can conclude that even if EWPT or QCDPT generated PMFs 
have a small amplitude at large scales, if 1 — 10% of the radiation energy density is in the 
form of a PMF, the observable CMB anisotropy consequences will be similar to those of 
an inflation-generated PMF with 10~ 9 — 10~ 10 G now at 1 Mpc. 

Another interesting effect is relic gravitational wave (GW) generation [ [H [9], [M] . The 
amplitude of these GWs is determined by the amount of MF energy density present at 
the PT, and it is not influenced by further MF damping. Direct cascade MHD turbulence 
generates GW with amplitude 

M/)-. 2 x 10 -W^^)(^) 1/3 (£) 1/2 ^ (/) , (3) 

where / is the frequency now, f H = Xjj 1 ~ 1.6 x 10" 5 Hz (^/100) 1 / 6 (T*/100 GeV) is the 
Hubble frequency now, and S(f) is determined by the MF statistical properties, 



S(f) = C 2 M I ^expf erfcf --rP^\ , (4) 



Re3/4 dx f /V 

where erfc(x) is the complementary error function defined as erfc(x) = 1 — erf(x), where 
erf(x) = f£ dyexp(— y 2 ) is the error function. As expected, the integral in Eq. (j4J) is 
dominated by the large scale (x ~ 1) contribution. 

The amplitude and the energy density of the GW are related through, 

h c (f) = 1.26 x 10- 18 (*j \ [hi n ow (/)] V2 , (5) 

where ficw(f) is the GW spectral energy density parameter. Using Eqs. (jSJ) and (jSJ), we 
have 

G w(/)^ = 2.3 x 10- 4 v 2 Al 5 (^) V3 3 S(f). (6) 

Integrating fi GW (/) over frequency, it can be seen that the efficiency of GW production is 
low, oc f^7 2 , with peak GW frequency for the EWPT being / pea k — I^va^h 1 [S5] (and 
an additional peak at A^ 1 for helical MHD turbulence [[36]), but the signal is potentially 
observable by LISA [ [M]. For lower frequencies / ^gw oc f 3 , and for high 

frequencies / ^> fn there is exponential damping. The peak amplitude of ficw(/peak) = 
2.3 x 10 _4 f^7 2 (^/100) -1 / 3 and is independent of T*. The analysis above shows that the 
main contribution to the GW background comes from the EWPT and we can ignore the 
expansion of the Universe when studying GW generation (even for a helical PMF). 
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4. Conclusion 

We have constrained a causally generated PMF produced prior to BBN, by using the 
BBN bound on the relativistic energy density during nucleosynthesis. We can also con- 
strain a PMF generated after nucleosynthesis, but still during the radiation-dominated 
epoch, by requiring that the PMF energy density not be the dominant component. Figure 
1 shows the PMF limits without accounting for damping or decay of the PMF. This is 
because PMF large-scale decay is very model dependent, and to study this case requires 
accounting for the initial conditions. 

We have also argued that directly using the smoothed PMF, B\, might cause confusion. 
We instead propose using the PMF energy density when deriving constraints from PMF 
cosmo logical signatures. In particular, using £ m (i) oc t n « (where the correlation length 
is associated with the largest size MF eddy A = 2rr/k ) together with £ M oc t~ nE , in 
the framework of £\f oc Bl(\k ) a+1 , implies that B x (t) oc £((<*+iK-^s)/2 ; leading to B x 
increasing if the smoothing scale A is fixed, while the real MF energy density is decreasing. 
It is obvious that B\ increasing in time is not a physical effect. Another advantage of 
using the PMF energy density is the independence of the BBN limit on the energy scale 
when the generation occurs (T*), even though all characteristic length scales are strongly 
T* dependent. 

Free turbulence decay on large scales (even without knowing the exact scaling law) 
allows us to conclude that after the PT ends turbulence has been largely damped, and 
there is no available source to produce the gravitational radiation at the same level as that 
from the PT itself. This allows us to neglect cosmological expansion when considering the 
GW generation process. Accounting for this, the direct detection of relic GWs will allow 
us to study the PT MHD turbulence picture, if enough (1-10%) of the thermal energy 
during the PT is present in the form of MF energy density. 

We appreciate helpful comments from A. Brandenburg and K. Jedamzik and acknowl- 
edge useful discussions with C. Caprini, R. Durrer, S. Huber, L. Kisslinger, A. Kosowsky, 
T. Stevens, K. Subramanian, and T. Vachaspati. We acknowledge partial support from 
GNSF grants Pres-07/153 and ST08/4-442 and DOE grant DE-FG03-99EP41043. T.K. 
acknowledges the ICTP associate membership program, and NORDITA for hospitality 
during the Electroweak Phase Transitions workshop. 
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